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Abstract. 

We consider the covariant Lyapunov vectors (CLV) of a high-dimensional 
Hamiltonian flow in the case of long range potential, namely the Hamiltonian 
Mean Field (HMF) problem, by studying the behavior of the Lyapunov spectra 
and the Oseledets' splitting principal angles (the mutual orientation between 
stable and unstable subspaces) when a phase transition takes place. Motivated by 
several results connecting the dynamical properties of systems to their Lyapunov 
exponents and vectors, we first find confirmation of an explicit sensitivity of such 
quantities to the transition: our main finding is that the mutual orientation 
between stable and unstable subspaces fluctuates in time around specific constant 
values. Moreover, the fluctuations statistics are different above and below the 
critical threshold, and thus intimately connected to the transition. Pictorially, the 
evolution of the stability subspaces along a reference orbits turns out to happen 
only through rigid global rotations plus weak internal deformations between the 
two subspaces; we are not aware of any theoretical explanation for such a non- 
trivial setting but, through analogies with integrable systems, we provide some 
little insight on 'how much' the system resembles an integrable one. 
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1. Introduction 

Through this paper we focus on a paradigmatic example of long-range interacting 
systems, the so called Hamiltonian Mean Field (HMF) model, first introduced by 
Antoni and Ruffo in 1995 [5] in order to model the behavior of galaxies or plasmas. The 
HMF is a toy-model that was specifically designed in order to be easily implemcntablc 
and theoretically studied, thus making it the ideal candidate to investigate the 
properties of long-range systems in the thermodynamic limit. Since its introduction it 
has been extensively studied |13) , and has become a paradigmatic model for the study 
of out-of-equilibrium quasi-stationary states (QSS). 

The HMF shows an interesting dynamic and thermodynamic behavior: when 
performing numerical simulations, starting out-of-equilibrium, the system is usually 
trapped in long-lived QSSs, before relaxing to the Boltzmann equilibrium solution. 
The QSSs exhibit a breakdown of ergodicity, a distribution of momenta which is 
non-Maxwellian and anomalous diffusion [9l [101 Qj]. The lifetime of these states 
diverges with the number of degrees of freedom according to a power law, so in the 
thermodynamic limit they become proper equilibrium states. The simplicity of the 
model allows one to develop a theoretical approach that explains the presence of QSSs 
as equilibrium solutions of the Vlasov equation, as pioneered in [T^ . However we must 
note that any finite-size HMF system will eventually thermalize to the Boltzmann 
equilibrium in a finite time for any choice of initial conditions. 

The HMF was also proposed to interpret the behavior of wave-particle interacting 
systems: under specific conditions its nonequilibrium properties reproduce the 
dynamical aspects of the Collective Atomic Recoil Laser (CARL) [6]. 




Figure 1. Pictorial presentation of the HMF system, where particles are 
identified by their angular position qj over a circle. The interaction is only 
sensitive to the angular distance A.qj^ between each pair of particles. 
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2. Framework 



2.1. The Hamiltonian Mean Field 

The HMF model describes TV identical particles interacting through cosine potentials: 



JV 
i=l 



N 

\p1 + m X^ 1 - cos (* ~ &)) 

3=1 



(1) 



The canonical coordinate qi identifies the position of particle i over the unitary circle 
[0, 2tt[, with periodic boundary conditions, while pi is the conjugated momentum; since 
the interaction depends on Aqjk = qj — qu only (see figure [T]) the system is invariant 
under global rotations around the unit circle and thus the total angular momentum 
p = i ^2jLiPj is constant (as H) along any trajectory. Notice that, moreover: 

• particles have the same mass, which we assume equal to one for simplicity; 

• particles can cross each other, or, equivalently, do experience perfectly elastic 
collisions, in agreement with the hypothesis that masses are identical. 

Depending on the value of the coupling constant C, the interaction can be either 
attractive (C > 0) or repulsive (C < 0); both cases were studied in the past [l3 | IT4 ] : 
in our computations we will consider the case (7 = 1. The HMF potential can be 
also seen as the first harmonic in the expansion of the one-dimensional self-gravitating 
potential with periodic boundary conditions [151 116j . 
Defining ' = ^ and V q = ^ , the canonical equations of motion read: 

q = P , 

JV 

P = - V q V = q q 3 = %Yl sin (* _ ( 2 ) 

fc=l 

or, denoting x = [q, p] and the velocity field v(x) = [p, — V q V], dynamics is given by: 
x = v(x) . (3) 

In this case it is then possible to formally write the Hamiltonian flow ft , such that: 

ft(x) = vof t (x) , f t of T = f t+T , f = id (4) 

for every point x in the phase space fi, an orbit is the parametrized curve x t = f{(x); 
by definition, velocity field is everywhere tangent to such curve, being the natural 
direction along which perturbations do not propagate in time. 

2.2. Stability 

The stability of flow Q can be investigated by considering infinitesimal perturbations: 
<Jx t =F t (x)$x , (5) 
where the Jacobian matrices cocyle is defined as: 

F *W^^ f t | x ; (6) 

obviously Fo(x) = 1 for all x, and the cocycle composition property is sathisfied: 

F t (x r )F T (x) =F t+T (x) . (7) 
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By the latter and from the symplectic structure of Hamilton evolution we get: 



Ft(x)" 



F-t(xj) = YF((x) T Y 5 



I 

1 O 



(8) 



where Y is a representation of the fundamental symplectic form. By using and 
© we get: 

F t (x) = K(x t )F t (x) , (9) 
with K the stability matrix, defined as the Jacobian matrix of velocity field v: 

K ( x )^ JHx . K(xt) = K t . (10) 
By deriving (j&J) w.r.t. time t, we obtain the equations of motion for the perturbations: 

5x 4 =K(x f )<5x t ; (11) 

tangent to fi in x t by definition and evolving exactly as each column of matrix F t . By 
its structure, both as a matrix and as a function of phase space, K induces the formal 
solutions to © as time-ordered exponential integrals, defined by: 

t/r 



F t (x) = T-exp 



/ drK(x T ) = lim ff [l + rK fer 
In t->0 j - j - 



(12) 



with the proper matrix-product ordering; notice that, in general, the exponent is not a 
time-average of the stability matrix unless the latter commutes with itself at different 
points of the chosen orbit; in such a case we indeed have: 

[K t ,K T ]=0 => Ft =c'< K ^ x ) (13) 

with [•, •] the usual matrix commutator and (-)t(x) the time-average for r G [0,i\. 
Such simplification always occurs in the integrable case of linear velocity fields and, in 
general, when considering the linearized dinamics of a fixed point; in both cases the 
stability matrix is constant along the chosen reference orbit (that is, any orbit for the 
constant system and the fixed point itself ) and thus sathisfies (|T3j) : 



K(x t ) = K(x) = K =>■ F f = c 



51 17\ 



t n K r ' 



(14) 



n=0 



In this case we can compactly express the spectral problem: we have to find a pair of 
matrices A , E , with A diagonal, such that : 

O I 

-H O 



KE = EA 



K 



H 



d 2 v 

dqi- l dq' y 



(15) 



with O , I and H = H T respectively the null, the identity and the Hessian N x N 
matrix of potential V and where we have used ([2]), ([3]) and (IT0| . The solution of (TT5"]) 
can then be expressed in terms of H as follows: 

A A 
Ax Ax 



E 



, A = 


X 


o 




o 


-x _ 



(16) 



with x = V— S diagonal and A, £ respectively the eigen-vector/values of matrix H: 
HA = AS , 
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Thus the eigenvalues Xk = V^^k are purely real/imaginary precisely when the real 
eigenvalues , correspondent the principal curvatures of the potential energy surface, 
are respectively negative/positive; this can be stated by rewriting (fT4|) as: 

F t = c tK = Ec IA E- 1 (18) 

showing that E is an eigenbasis for F t at all times, with diagonal action e tXk , for each 
k, corresponding to real exponential contractions/dilations (negative curvatures) or to 
complex rotations (positive curvatures). This characterizes dynamics (1111) completely: 
the subset {e^ } of real vectors of basis E corresponding to negative curvatures £& < 
are exactly the physical directions along which the system equally contracts and 
expands, in agreement with its area-preserving nature; indeed, by choosing: 

6x = e± => dx t = e ±tXfc e± (19) 

one can see that the set of Lyapunov exponents A arise [2] : 

hm A(x,c5x) = flog&ji =± X k =±A« • (20) 

t^oo 1 ||dx|| 

Nonzero Lyapunov exponents are associated with the subset of real eigenvalues of K; 
this yields the Lyapunov spectrum (see figure [2] for an example); the positive curvature 
directions, associated to imaginary eigenvalues ±i|Xj| instead, extend in the complex 
domain and yield the principal axes for ellipses in each of the planes (qj,Pj) invariant 
under the associated rotations of frequency \xj\- 
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3. Oseledets' Theorem 
3.1. Transformed Cocycles 

In the generic non-constant nor fixed-point cases, however, the eigenbases of K(x) 
evaluated in points along orbits do not diagonalize cocycle F t , as shown by the relation: 

ft 



F t (x) =E(x t )T-exp 



J dr (a - e^e) ( Xt ; 



E(x)- 



(21) 



and (f2Tj) reduces to (|18|) only when E = O, i.e. the eigenbasis is constant along orbits: 

E(x) eA[ E o f t ] t=Q (x) ee 2? X E - O ; (22) 

here £> x is the material derivative for phase-space dependent observables, or, the 
derivative along the flow: 

27V JV 

^£^)at--E(^4-C4)- (23) 
fc=i fc=i 

This fact naturally induces to consider the problem of finding a phase-space dependent 
basis U(x) generating a transformed cocycle F t (x) : 



F t (x) = U(x t )F t (x)U(x)" 
and a corresponding transformed stability matrix K(x) 



F t (x) =T-exp 



drK(x r ; 



(24) 



(25) 



constrained to assume some prescribed shape; for example, to obtain a characterization 
like in (TTg)) for generic systems, one would ask for a diagonal transformed cocycle; 
combining equations (f2~Tj) and (|2"4"|) for a given phase-space dependent basis U yields: 



F t (x) =T-exp 



I dT ( u ~'( 



KU-U (x T ) 



(26) 



and through (|25j) . (|26j) . once given a structure for matrix K, we obtain: 

KU - £> X U = UK , (27) 

that is, a PDE for U(x) defined over some (possibly dense) subset of phase-space; 
notice that equation (|2"T1) is completely local and time-independent, can have at most 
2N linearly independent solutions and, when asking for diagonal transformed cocycles 
and restricted to linear velocities or fixed-points, reduces to the eigen-probem (fTS")) . 

3.2. Covariant Lyapunov Vectors 

The basis U producing a diagonal transformed cocycle from F t is called its Covariant 
Lyapunov vectors basis (CLV), or Oseledets' splitting, and is associated to the diagonal 
transformed stability matrix K = A so that equation (|27|) takes the form: 

KU - 2? X U = UA . (28) 

and A(x) is the set of local exponents whose time-average determines the new cocycle: 



Ft(x) 



2 tA '( x > , A t (x) = <A) t (x) = \ [ drA(x T 

Jo 



)• 



(29) 
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In the framework of nonvmiformly hyperbolic systems, the Oseledets' theorem [3] then 
states the existence of such precise decomposition of cocycle F t into a phase-space 
dependent 'basis part' U and a time-dependent 'exponential part' A t , such that: 

F t (x)=U(x t )e tA *( x )U(x)- 1 , (30) 

where the exponent matrix satisfies the infinite time limit (Kingman's theorem, [5]): 

lirxi A t (x) = A{x} = (AMx) . (31) 

t— >±oo 

with {x} = {fi(x)} t6 R the orbit of x, along which A{x} is constant, or, flow-invariant. 
Thus the tangent dynamics exactly maps each vector field (x) (a column of U) on 
itself evaluated at a new point, u^(x t ), expanding it by the same asymptotic exponent 
A/, in both temporal limits [2]; the quantity is called the Lyapunov exponent 
associated to the fc-th covariant Lyapunov vector (CLV), while the sets of all exponents 
and vectors go respectively under the names of Lyapunov spectrum A and Oseledets ' 
splitting U of cocycle F, usually ordered by the rule: 4=> Aj > A*, (as in figure[2]). 
Throughout this paper, positive local exponents will be generically referred to as 
Ai, with the label 1 standing for 'one-step' due to the fact that, numerically, such 
quantities are calculated at every step of the CLV algorithm [24]; in agreement with 
the second equation in (|2"5|) . the FTLE A calculated over a time t will be labeled (Ai) t . 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

j/2N 



Figure 2. HMF, JV = 256, t = 10 7 : Lyapunov spectrum obtained by Benettin's 
algorithm [7] for the energy density value e = 0.65; notice that the spectrum is 
symmetric w.r.t. the normalized index j = N, as symplectic property J5J implies. 



3.3. Remarks 

3.3.1. Ordering In the generic Hamiltonian case, because of the symplectic structure, 
areas in each (gfc,pfc)-plane are conserved by the flow and the exponents are usually 
constrained to obey the symmetry: h.2N-k+i = — A& = i..jv; notice that this condition, 
as well as the exponents ordering named above, is not fulfilled by the present 
calculations since, by the convention used here, the matrix of local exponents reads: 



A = 



A O 
O A 



(32) 
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with A the diagonal matrix containing the positive exponents only; imposing the 
decreasing ordering on the latter induces an increasing ordering on the negative block, 
but this is completely unharmful, being equivalent to a global basis permutation and 
thus factorizable over each considered matrix. In the numerics we apply decreasing 
ordering to the spectrum, obtaining match with the symplectic condition (figure [2]). 

3.3.2. Numerical Approximations The search for upper/lower triangular cocycles 
(implying upper/lower matrices K) yields orthonormal phase-space dependent bases. 
The latter are called Oseledets' orthonormal bases and can be usually approximated 
by Benettin's algorithm [7], which by-produces an approximation of the Lyapunov 
spectrum also; such bases are the same used to perform Ginelli's algorithm [21] to 
approximate the whole splitting. 

3.3.3. Tangencies It should be noticed that the possible existence of solutions U of 
equation ([27)) crucially depends on whether the condition det U ^ is fulfilled or not, 
i.e. whether the corresponding column vectors do form a basis; for chaotic Hamiltonian 
systems this is supposed to hold in the whole phase-space exept for subsets of zero 
Lebesgue measure, where two or more CLV can coincide, making the matrix U locally 
non-invertible. Yet, it is not clear how the effect of such pathological sets, usually 
called tangencies, can affect the dynamics of systems with a large numbers of degrees 
of freedom; the main problems are connected to the smoothness properties of the flow, 
that naturally induce the splitting to arbitrarily shrink in points closer and closer to 
the tangency. In such regions outside tangencies the splitting is still invertible but in 
a numerically weak sense, i.e. all the covariant directions are transversal but some of 
the principal angles between them may be very close to zero (see paragraph I4.2[) . 
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— Theoretical prediction 



» N=512 
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Figure 3. Second order phase transition between homogeneous and magnetized 
states for the HMF system at equilibrium, at a critical energy per particle 
e c = 0.75. The predicted theoretical curve, displayed by the red thick line, is 
obtained from the minimization of the free energy in the canonical ensemble. It 
is consistent with the numerical microcanonical simulations carried out for finite 
N (the positive magnetization observed in the homogeneous phase is of the order 
of the finite-size fluctuations ~ ~7s)- 



4. Phase transition 

In the following we will focus on the thcrmodynamical equilibrium description of 
the HMF model, by studying finite-size systems and avoiding any potential QSS. 
The simplicity of the HMF allows to obtain analytical solutions for its statistical 
equilibrium properties, both in the canonical [5] and microcanonical ensemble [13 , 
independently and without relying on ensemble equivalence. 

When considering an attractive interaction, the HMF displays a second order 
ferromagnetic phase transition at equilibrium, which separates homogeneous from 
inhomogeneous states. In order to characterize this transition we consider the order 
parameter /x, which is usually referred as the magnetization, in analogy with spin 
systems, defined as the modulus of the vector M with axial components M x M y as: 



By defining e — E/N as the energy density of the system, the critical value that 
identify the transition is e c = 0.75 (see figure [3|), above which magnetization is zero. 
To study the transition we consider a sampling of 13 values for the energy density: 

• below threshold, e from 0.5 to 0.7 with a step Ae = 0.05 (blue/dashed curves); 

• near threshold, e = 0.73 (green/dashed curves), 0.77 (purple/continuous curves); 

• on threshold, e = 0.75 (black/ continuous thick curve); 

• above threshold, s from 0.8 to 1 with a step Ae = 0.05 (red/continuous curves); 

the described association between colors/line- types and energy densities is kept 
constant throughout the paper, as well as the choice of the system parameters: the 
number of particles, N = 256, and the orbits duration, t = 10 7 , in natural time units. 




(33) 
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Figure 4. HMF, TV = 256, t = 10 7 : total probabilities of all TV — 2 positive local 
exponents Ai for the 13 energy densities described in chapter [3] two exponents 
are always due to invariance of the Hamiltonian H and total momentum p. The 
distributions have compact support in all cases, but strong modifications occur 
along the transition: below trcshold the range of local exponents is broader (being 
higher the largest LE, see figures in IBI and 111 I for j = 1), giving rise to a plateau 
with several maxima, disappearing as the critical energy is approached; above the 
transition probability accumulates on the only maximum remaining, located in 
Ai = 0. This is already a rough argument supporting the fact that the disordered 
phase is far less unstable than the orderd one 1171 , although arbitrarily small 
FTLE never disappear from the latter (see again figures [51 1111 for j = TV). Notice 
that negative Ai values in these pictures are due to fluctuating violations of the 
symplectic condition J5J; such fluctuations, although very small and decreasing as 
1/t, are stronger and last longer when the stability is weaker, i.e. above threshold. 

4--1- Lyapunov spectrum 

As a first test to probe the sensitivity of tangent dynamics on the phase transition, we 
consider the distribution -P(Ai) for the set of first N — 2 local exponents X%, i.e. the 
nonzero ones, discarding for now the information associated to the spectrum label j. 
By inspecting figure[4]we notice that some qualitative features change between the two 
phases: while in the homogeneous phase the distribution looks like an antisymmetric 
bell around Ai = 0, an instability plateau appears below the critical threshold; in all 
cases the distribution have compact support, by definition, over Ai positive values. 
It was as already noticed [17] . from numerical analysis, that the first local exponent 
exhibits anomalous fluctuations when below threshold, increases with energy up to e c , 
while for e > e c it goes to zero. The time averages of Ai for each value of the FTLE 
spectrum index j are then organized in figure [5j showing that the ordered phase has 
larger exponents than the disordered one, although arbitrarily small values persist 
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Figure 5. HMF. TV = 256, t = 10 7 , positive Lyapunov spectra for the considered 
13 as described in chapter [4] a steep decay in the instability is apparent exactly 
above the transition e = 0.75, with a decrease of the whole spectrum and 
important attenuations on the larger exponents. Arbitrary small exponents persist 
at all phases and are tipically responsible for anomalous statistical features of the 
dynamics, as shown by Ai variances in the next figure. 




Figure 6. HMF, N = 256, t = 10 7 : time variances of the positive local 
exponent Ai for 13 values of the energy density e as described in chapter [4] 
such quantities measure the broadness of fluctuations (decaying as 1/t) of each 
FTLE (see also figure [Til l. The transition changes the exponents statistics from a 
more uniform situation, where the fluctuations over all exponents are comparable 
(magnetized phase, blue/dashed curves), to a situation in which the small and 
intermediate exponents fluctuates much less than the large ones (disordered phase, 
red /continuous curves) . 

at all phases. Fluctuations in time of the FTLE spectra are also sensitive to the 
transition, as shown by the variances o"(Ai) in figure in the disordered phase, large 
exponents have much stronger fluctuations than the intermediate and small ones, while 
in the magnetized phase the exponents tend to fluctuate more uniformly. Finite size 
effects are also important: in the limit of infinite size N — > oo the largest Lyapunov 
exponent (LLE) drops to zero with a power law scaling N~s [TS]. An analytical 
estimate of the scaling of the LLE was proposed [23] j which uses a geometrical method 
based on microcanonical averages [TS], and the same law was also obtained from a 
random matrix approach |21j ; the model was shown to move towards integrability for 
N large also in [35] . A study of the scaling properties of the Lyapunov spectrum and 
of the Kolmogorov-Sinai entropy was performed in [20) . 
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Figure 7. HMF, N = 256, t = 10 7 : Principal angles for 13 values of the 
energy density e as described in chapter[4] showing a qualitative change across the 
transition, with persistency of arbitrarily small angles above and below e c = 0.75. 

4-2. Principal Angles & Hamiltonian Stability 

We now proceed to illustrate our main numerical findings, by probing geometric 
features of tangent dynamics. As already mentioned, when some of the covariant 
vectors take very similar directions it is possible that, nearby, there exist points in 
which the splitting is not invertible, i.e. tangencies appear. In general, given two 
sets both spanned by n < 2N unit vectors in M. 2N , the reciprocal orientation between 
the two is defined by the so called principal angles (PA), defined recursively, from the 
smallest to the largest, as the smaller possible angles between the two subspaces [4]. 
More explicitly, writing the two sets of vectors as 27V x n matrices U and S, such 
angles then correspond to the singular values of the matrix product: 

U T S = AZB T , Z = diag(oos(0 J -))i=i...» (34) 

with the singular values decomposition (SVD) of product U T S on the r.h.s. of the first 
equation; the two basis matrices A, B are orthogonal and, along with the diagonal 
matrix Z of real eigenvalues, they diagonalize the two symmetric matrices: 

U T SS T U = AZ 2 A T , S T UU T S = BZ 2 B T (35) 

In the present work, in order to probe the hyperbolicity of the system, we compute the 
PA between the stable and unstable CLV, i.e. the pair of 2N x N matrices composing 
by blocks the CLV matrix, U = [U+, U_] . In figure [7] the spectra of principal angles 
are shown for the 13 values of the energy density already described in this chapter; 
we notice a strong qualitative change in the angles profiles, with a flip of the sign 
of the concavity as the main feature. On the other hand, both above and below 
the transition there are arbitrarily small angles, meaning that the system remains 
nonuniformly hyperbolic for all energy densities; in other words, tangencies or almost- 
tangencies are always present between the stable and unstable subspaces, in agreement 
with the fact that for all energy values there are arbitrarily small FTLE (see figure 
Notice also the extremely low amplitude of fluctuations (see figure Q2] for an overlook); 
this feature is strong in the ordered phase. This is a very subtle property: apparently, 
tangent subspaces respectively spanned by the stable/unstable CLV correspond to 
structures whose reciprocal orientation is very rigid. Practically speaking, this means 
that the splitting changes along a reference orbit through arbitrary global rotations 
(i.e. all the vectors rotate through the same matrix, leaving the splitting undeformed) 
plus very weak internal deformations. 



Oseledets ' Splitting in Large Hamiltonian Systems: the HMF Phase Transition 13 




Figure 8. HMF, N = 256, t = 10 7 : time variances of the principal angles for 
13 energy densities as described in chapter [4] and for all the angles {dj}j=i...Ni 
notice the huge qualitative and quantitave variation across the transition: below 
the threshold, in the disordered phase, all the fluctuations are very low, while 
exactly across the transition (e = 0.75, black/continuous thick curve), and above, 
the spectrum fluctuates by one order of magnitude more; interestingly, in the 
disordered phase (e < 0.75) there can be identified a subset of the spectrum, 
j € that shrinks to zero while moving toward the transition and in which 

the angles (from the largest to 0j*) have very low fluctuations; although surely 
connected to the geometrical aspects of the transition, this last fact remains subtle 
to explain due to its statistical nature. 



4-3. Principal Angles for Fixed Points 



The small fluctuations of principal angles suggest to explicitly calculate them for fixed 
points, with K constant and U = E as in (fT6| and (fT?]) ; by taking into account that 



each column of E is normalized to 1 by the factor 1 + 5L , decomposition (|34j reads 



E T E 



W,- 



U T IF 



U T U_ 





I 


w 




w 


I 



1 



i-x 



1 



(36) 



(37) 



with W diagonal, \j a rea l eigenvalue of K and Ej the corresponding negative 
principal curvature/eigenvalue of the Hessian matrix H of potential V; since the 
singular decomposition of each (in this particular case, diagonal) block matrix yields 
the identity matrix as orthonormal bases, matrix W contains itself the singular values 
and thus the cosines of pricipal angles between stable U_ and unstable U + subspaces: 



cos( 



(38) 



yielding a relation between principal angles and local exponents, that in this case 
correspond to the exact Lyapunov spectrum Aj. 

Such relation between principal angles and Lyapunov spectrum is a very interesting 
property that, in principle, can be tested over any Hamiltonian system once the spectra 
of exponents and angles are known; since the relation holds exactly for potentials that 
are quadratic functions of the configurational coordinates q, here we quantify how 
much the system deviates from a quadratic one by the difference: 

Ri = cos( ) - 1 - ; (39) 

1 u 3 " 1 + (A{) 2 V ; 
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Figure 9. HMF, N = 256, t = 10 7 : the difference R in equation g9]l for the 
13 values of the energy density as described in chapter [4] the FTLE spectrum, 
obtained by time averages (X\)t, is treated as the set of eigenvalues of some 
Hessian matrix H that is constant in phase space, i.e. corresponds to some 
potential V that is quadratic in the q coordinates. When the system is truly 
quadratic, R = for all the spectrum; here we find that, above threshold (e > e c , 
purple and red continuous curves) large subsets of the FTLE and principal angles 
spectra behave as in a quadratic system. 

obtained by the time-averages of both the angles spectrum and the local exponents 
spectrum; figure [9] shows that the system never coincides with a purely quadratic one 
but, once above threshold e c = 0.75, becomes nearer and nearer to it, with larger 
and larger portions of the spectra having R falling near zero. This can be seen by 
plotting the RMS of R (over j) versus the order parameter, figure [TUl fitting well a 
linear profile with slope ~ 1 that, interestingly, localizes the RMS vanishing point at 
[i = 0. These facts suggest that, in the ordered phase, the tangent structure of the 
system is not uniform across phase space although the reciprocal orientation between 
stable/unstable subspaces is rigid; on the contrary the disordered phase, characterized 
by larger internal fluctuations, has a large subset of tangent subspaces that behave as 
in a uniform, quadratic system. 
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Figure 10. HMF, N = 256, t = 10 7 : scatter plot of the the root mean square of 
R over j l|39|l VS the order parameter /i, for the 13 values of the energy density 
as described in chapter [4] with the correspondence: below threshold (dots), on 
threshold (cross) and above threshold (circles); we find good agreement with a 
linear fit (red line) with slope ~ 1, that gives a vanishing RMS point at fi = 0. 
Notice that in all the performed simulations the resulting RMS is nonzero, in 
agreement with the fact that the system under study never generates a true 
quadratic effective potential but, as shown here, deviates from it by a quantity 
highly correlated with the order parameter of the system. 



5. Conclusions 

We considered the Hamiltonian Mean Field problem for a range of energy density 
e from below to above the analytic critical value e c = 0.75, calculating for each 
case both the Lyapunov exponents and the Oseledets' splitting principal angles 
spectra; the sensitivity of the whole Lyapunov spectrum to the phase transition is 
confirmed, with strong correlations with the behavior of magnetization // (the order 
parameter of the system). Our main result is that a dynamical analysis in terms 
of the Oseledets' splitting angles turns out to provide a clearer dynamical signature 
of the phase transition, with respect to more conventional indicators based on the 
Lyapunov spectrum only. Through extensive numerical investigations, we find that 
stable/unstable subspaces in the /j, ^ phase obey to a nontrivial constraint for which 
the two tangent subsets (evaluated pointwise along some reference orbit) mantain a 
rigid reciprocal orientation that can only fluctuate very weakly in time. 
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Figure 11. HMF, N = 256, t = 10 7 : time -average (black curves) and time- 
averageivariance (red curves) of local exponent Ai spectra for the considered 
energy density values: e = 0.50, 0.55, 0.60, 0.65, 0.70 (resp. from upper panel 
in left column), e = 0.73, 0.75, 0.77 (resp. from upper panel in central column) 
and e = 0.80, 0.85, 0.90, 0.95, 1.00, (resp. from upper panel in right column), 
i.e. the same graphs of figure [5] with superimposed the fluctuations obtained 
by the variances in {6}; the phase transition (e c = 0.75, central panel) induces 
modifications both in the maximal value of the spectra (in j = 1 or j/2N — > 0, 
by definition) and in the behavior of the corresponding variances. Above the 
transition (all right column and lower panel in central column) fluctuations are 
overall wider than below, while the maximal instability is weaker. 
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Figure 12. HMF: time-average (black curves) and time-variance (red curves) of 
principal angle spectra for the considered energy density values: e = 0.50, 0.55, 
0.60, 0.65, 0.70 (resp. from upper panel in left column), e = 0.73, 0.75, 0.77 
(resp. from upper panel in central column) and e = 0.80, 0.85, 0.90, 0.95, 1.00, 
(resp. from upper panel in right column); the transition is apparent through the 
variances of the spectra: indeed, for energies below threshold (e < e c = 0.75, 
all left column and upper panel in central column) the fluctuations around the 
average are very small. Above threshold, instead, fluctuations rise significantly, 
while the concavity of the average spectra flips its sign; notice that in all dynamical 
regimes the largest /smallest angles share small fluctuations, indicating that the 
upper /lower bounds for the splitting amplitude are practically constant and 
the fluctuations in the geometry of stable/unstable subspaces affects mainly the 
intermediate angles. 
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